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Abstract. We examine the nature of the cosmological big-bounce (BB) transition within the loop 
geometry underlying Loop Quantum Cosmology (LQC) at classical and quantum levels. Our canon- 
ical quantization method is an alternative to the standard LQC. Our method opens the door for an- 
alyzes of spectra of physical observables like energy density and volume operator. We find that one 
cannot determine the energy scale specific to BB by making use of the loop geometry without an 
extra input from observational cosmology. 
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INTRODUCTION 

Observational cosmology strongly suggests that the universe emerged from a state 
with extremely high energy densities of physical fields. Mathematical cosmology gives 
numerous examples of solutions with an initial Big-Bang singularity. This is why it is 
commonly believed that the cosmological singularity is a real feature of the Universe. 

It seems that one can treat the problem of the singularity by making use of loop 
geometry: classical Big-Bang turns into quantum Big-Bounce transition. Thus, one may 
try to answer some interesting questions: What is the energy scale of the unification 
of general relativity (GR) with quantum physics? What is the structure of spacetime at 
semi-classical level (simply connected, foamy or discrete)? What is the origin of matter, 
non-gravitational fields and spacetime? What was before the Big-Bounce? 

We try to address the first two questions in the context of the fiat FRW universe with 
massless scalar field. 

There are two methods of quantization of the cosmology models based on loop ge- 
ometry: (i) standard loop quantum cosmology (LQC), based on the rule: 'first quantize, 
then impose constraintsrL2j, and (ii) non-standard LQC using the rule: 'first solve con- 
straints, then quantize In the former case, one believes that the classical Big-Bang 
is replaced by quantum Big-Bounce due to strong quantum effects at the Planck scale 
[H^laB]- In the latter case, one shows that it is the modification of GR by loop geometry 
which is responsible for the resolution of the singularity and the energy density of matter 
at the Big-Bounce is unknown Jsj]. 



CLASSICAL LEVEL 

The gravitational part of the GR Hamiltonian reads [H] 

Hg:= f d^xiN'Oi+N^Ca+NC), (1) 

where E, space-like part of spacetime R x L; {N',N",N), Lagrange multipliers; 
(C/,Ca,C) are Gauss, diffeomorphism and scalar constraints; {a,b = 1,2,3), spatial 
indices; (z, j,^ = 1,2,3) internal SU{2) indices. Constrains must satisfy the specific 
algebra. For flat FRW universe with massless scalar field (in some special gauge) the 
Hamiltonian ([U) turns into [H] 

Hg = -r-2 f d\ Ne-'eijkE'^^E^^Fi^^ , (2) 

where y, Barbero-Immirzi parameter; ^ C E, elementary cell; A^, lapse function; e,jyt, 
alternating tensor; Ef, density weighted triad; F^^^ = daA'l^ — d^A^ + ef^jA'^Aj^, curvature 

of SU{2) connection A^; e := A/[detE|. 

Modification of GR by loop geometry means replacement (]^ of by its approxi- 
mation 

Fi(A)=^lim{-2r.(-ii^).'X".i}, (3) 

and we have 

F„i= lim Ff,(A), (4) 

A — > 

where Vq is the so-called fiducial volume. Holonomy of the connection around the square 

1/3 

loop with sides length ^Vq is found to be 

h\f=h\>'^hf\h^^Y'{hfr' (5) 

1 /3 

The holonomy along straight edge "e^da of length iiVq (in fundamental, 
j = 1/2, representation of SU(2) group) reads 

hl^\c) =exp(T^jUc) =cos(jUc/2) 1 -f2 sin(;Uc/2) T^, (6) 

where = —iOk/2 (o^ are the Pauli spin matrices). 

Making use of the so-called Thiemann identity leads finally to 

where 

fff > = -^^IiV£y'r.(/,W/,J«(*<*')->(/.i«)-'/.l«{(/.f V'.v}). (8) 



3 T 

and where V = \p\^ = a Vq is the volume of the elementary cell y. Variables c and p 
determine connections and triads E^:A'^ = "co^cVq and E^^ = "e^^/qJ,pVQ '^^'^ , 

where c = yaVj'''^ and \p\ = o^Vq^^, {c,p} = SnGy/^. 

Total Hamiltonian for FRW universe with a massless scalar field ^ reads ^ 

H = Hg + H^, (9) 

where Hg is defined by dV]) and //^ = p'^\p\^i /2, and where and p^j, are elementary 
variables satisfying {^.p^} — 1. The relation H ^0 defines the physical phase space. 

Making use of we calculate ([8]) and get the modified total Hamiltonian corre- 
sponding to ^ 

" /^-~8^^^'+2^' ^^^^ 



where 



/3:=^, v:=b|3/2 (11) 



\P 



are the canonical variables (of so-called improved scheme). It should be emphasized that 
(fTOl ) presents a purely classical modified Hamiltonian (with no quantum corrections)! 
The Poisson bracket for the canonical variables (/3, v, is defined to be 



d^dp,j, dp<j,d^' 



dpdv dvdp 

The dynamics of a canonical variable ^ is 

|:={^,//W}, ^g{/3,v,<^),p4, (13) 

where := d^/dT, and where T is an evolution parameter. Dynamics in physical phase 
space, '^phys^ is defined by solutions to (fT3l) satisfying the condition H^^^ ^ 0. Solutions 

of (fT3]) ignoring the constraint //^'^^ ^ are in kinematical phase space, ^j^j^J ■ 

If the Hamiltonian is a constraint which may be rewritten as a product of a simpler 

(X] 

constraint and a function on with no zeros, the original dynamics may be reduced 
to the dynamics with the simpler constraint Equation (flOl) can be rewritten as 

fjW ^^fjWfjiX) (14) 

where 

where = AnG/3. 

Hq^^ = if£ P(j, = = sin(A/3). In such case H^^^^ = 0, thus H^^^ equals identically zero 
so there is no dynamics. We exclude such pathological cases from further considerations 

and assume that Hq 7^ 0. 



(X) 

For functions / and g on ^^i^y^ we have 

^ = U,A^//i^)^W}=A^//^^)U,^W}, for ^W^O. (17) 

The relation 



8 dg NH^Q^\g,HW} {S:HW} 
may be rewritten as dH] 

df dg 



as H^^VO, (18) 



(19) 



Thus, in the relative dynamics one canonical variable may be used as an 'evolution 
parameter', and the dynamics is gauge independent (no dependance on A^). 

Since the relative dynamics is gauge independent, we choose N = 1/Hq 'to simplify 
calculations. Equations of motion read [i3|] 

^ = Krsgn{p0), P0=O, (20) 
/3 = -4;rGr^^^^^^^, v = 4;rG7vcos(A j8) sgn(sin(A j8)), (21) 

A 

H^^^ ^ 0. (22) 
Solution of relative dynamics is found to be IH] 

2v = Acosh(3Ks{(j)-(j)o)-\nAy (23) 

where s := sgn(;?,^), A := KYX\p(j,\. Solution for /3 may be determined from (|22l) rewrit- 
ten as H 

v| sin(A/3)| = K-yA Ip^I- (24) 
We can see that the variables v and j8 are functions of an evolution parameter 0. 

A function, : — > i?, is a Dirac observable (we choose N = 1/Hq as our 
method is gauge invariant) if 

{^,//(^)} = {ff^NH^^^H^^^} = {^,^(^)} = 0. (25) 



Thus ^ is solution to the equation dsj] 

^ T?r -^cos(Ap) — — ^— -=0. (26) 

A (9/3 ^ dv AkG (90 

Solutions to (|26l) are found to be lH] 

^i:=P^. ^2:=<^)-:r^arth(cos(A/3)), ^3:=,v^^^^^. (27) 



Observables satisfy the Lie algebra Jsj] 

{^2,^l} = l, {^i,^3} = 0, {^2, ^3} = 7'^. (28) 
Due to the constraint H^^^ =0,we have 

^3 = yK:^i. (29) 

Thus, in the physical phase space, '^p'^y^, we have only two observables which satisfy 
the algebra 

= 1. (30) 



Our kinematical phase space, ^^,„ , is four dimensional. In relative dynamics one 
variable is used to parameterize three others. Since the constraint relates two variables, 
we have only two independent variables. This is the reason we have only two observ- 
ables. 

We consider functions which can be expressed in terms of observables and an evo- 
lution parameter ^ so they are not observables. They do become observables for each 
fixed value of 0, since in such case they are only functions of observables. 

In what follows we consider the energy density of matter and the volume function. 
The energy density of matter field reads 

p(A,(|>) = ^^. (31) 
In terms of observables we have IH] 

p^ = ^i^ v=fC7A|^i|cosh(3fc((|)-^2))- (32) 
The density p takes its maximum at the minimum of v 



1 l2 



(33) 
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Let us apply ([33]) to the Planck scale: //>/ := ^fiG/c^ ~ 10 ^^m; ppi := c^/hG 
lO^'^GeV . Equation (l33l) in terms of suitable units rea 

3c2 1 

(34) 



Substituting A = Ipi into dH]) gives p,nax/ ppi - 2, 07; for pmax = ppi we get A ~ 1 , 44 //>/. 
Surprisingly, the classical expression (|34l) fits the Planck scale. 

The volume function reads: V(0) = v. It results from the constraint equation 

v|sin(A/3)| = K-7A|p,|,| (35) 



We use 7~ 0.24 determined in the black hole entropy calculations. 



that at the classical level the volume function is bounded from below [Ht] 

Vmm = KY^\P^\- (36) 

For A = the minimum of v vanishes! An interesting question is: Does the quantum 
volume operator have a non-zero minimum eigenvalue? 



QUANTUM LEVEL 

Our quantization method of constrained systems is simple enough to be fully controlled 
analytically. In the Schrodinger representation (since ^i, ^2 ^ ^) we have 

^2 ^ ^2 := ^ , ^1 := -ihd^: (37) 

where xeR, and where the carrier space is L^{R) ■ Thus, the representation of (|30l) reads 

[^2, ^1] = ih. (38) 

The representation is essentially self-adjoint on Cq{R) C L^(i? 
In the representation (l37l) the energy density operator reads 

P^P-= ^9 • (39) 

^ ^ 2(?C7A)2cosh2[3fc((^)-x)] 

Solution to the eigenvalue problem 

p\ir = p{xo)\{r, (40) 
for fixed value of ^, is found to be [0] 

xir{x) = d{x-xo), p{xo) = !..2 — u2 \a ^^^"^ 

2[KYkY cosh^ (3k-(0 -xq)) 

Energy density has the same functional form as the classical one! The range of spectrum 
is: (0, 2(j^^)- We can determine pmax if we know A. But, A is a free parameter of our 
method! Thus, finding the critical density of matter corresponding to the Big-Bounce is 
an open problem [8]. 

Since we have V = v, where 

V = KjX l^il cosh[3fc((^) - ^2)] = I'fyA ^1 cosh[3K((^) - ^2)]! =: |w|, (42) 

we first examine the spectrum of w. The quantum operator corresponding to w is defined 
as follows ^ 

w\\f:= T<Y^^ cosh[3K-(0 - (^2)] +cosh[3K-(0 - (^2) 
= +2cosh[3K:((^-x)]^-3K:sinh[3K:(<^-x)] ) (43) 



where \{f G D{w). The solution to the eigenvalue problem of the operator w reads jst] 

w\ifb = b^fb, beR, (44) 

./^exp(/^^ arctane3'^(^--^)) 

V/z,W := -^^ r ■ (45) 

cosh3 3fc(^ — x) 

For the volume operator, corresponding to V = v = we have 

y -^Vfb = \w\fb:=\b\fb. (46) 

Thus, the spectrum of the operator V seems to be continuous, and b — Qis the smallest 
eigenvalue of V. However, it turns out that the spectrum is discrete 

MINIMUM LENGTH PROBLEM 

Determination of A by standard LQC means [6]: 

• considering eigenvalue problem for the area operator, Ar = |;? | , in kinematical phase 

space of LQC: Ar|/i >= |/i >=: ar(/i) |/i >, where ar(/i) are continuous 

since HER; 

• making reference to discrete eigenvalues, {0, A, . . .}, of kinematical Ar of LQG, 
where A := 2^/3 Ttylp', 

• assuming: ar{X) = A, which leads to A = 3^/3/2. 

One assumes in the standard LQC that a surface cannot be squeezed to the zero value due 
to the existence in the universe of the minimal quantum of area. Physical justification 
for this assumption, offered by standard LQC, is doubtful because 

• Ar was examined so far only in kinematical Hilbert space of LQG, i.e. spectrum of 
Ar ignores the algebra of constraints of LQG so it has no physical meaning; 

• discrete spectrum of LQG was used to replace continuous spectrum of LQC (spec- 
tral discretization 'by hand'); 

• LQC is not a cosmological sector of LQG, but a quantization method inspired by 
LQG. Thus, LQC and LQG are models of different systems. 

The assumption of the standard LQC that low-lying eigenvalue of the area operator of 
LQG determines the basic parameter of LQC is ad /loc MIToll. 

CONCLUSIONS 

The modification of the classical Hamiltonian realized by making use of the loop ge- 
ometry, parameterized by A, turns Big-Bang into Big-Bounce. However, A is a free 
parameter to be determined. 



The determination of A is presently a great challenge. It seems that the FRW model is 
useless in this context; no eigenvalue of the volume operator is privileged. An extension 
of study to homogeneous (Bianchi I, ...) and isotropic (Lemaitre, ...) models cannot 
probably change the situation. It seems that it is the observational cosmology which may 
bring some resolution to the problem. So far, the cosmic photons reveal no dispersion 
[|ll|] up to the energy scale 5 x 10^^ GeV. One hopes that the primordial gravitational 
waves have made some imprints on the CMB spectrum so they may be used to determine 
the allowed range of values for the parameter A . 
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